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I/Tj ■ Abstract 

■ We study the capability of the international linear collider (ILC) to probe 

extra dimensions via the seesaw mechanism. In the scenario we study, heavy 
Kaluza-Klein neutrinos generate tiny neutrino masses and, at the same time, 
have sizable couplings to the standard-model particles. Consequently, a Kaluza- 
Klein tower of heavy neutrinos (N) can be produced and studied at the ILC 
through the process: e + e~ —> vN followed by N — > W£ decay. We show 
that the single lepton plus two-jets final states with large missing energy from 
this signal process will provide a good opportunity to measure the masses and 
cross sections of Kaluza-Klein neutrinos up to the third level. Furthermore, 
the neutrino oscillation parameters can be extracted from the flavor depen- 
dence of the lowest-mode signals, which give us information about the origin 
of low-energy neutrino masses. 



1 Introduction 



The detailed structure of lepton sector has been gradually revealed by the recent 
neutrino oscillation experiments [H |2]. The smallness of neutrino masses is one of 
the most important clues to find new physics beyond the standard model (SM). The 
seesaw mechanism naturally leads to small neutrino masses through heavy particles 
being coupled with ordinary neutrinos. In Type I seesaw scheme [3], the introduction 
of right-handed neutrinos implies intermediate mass scales to have light neutrino 
masses of order eV, and hence these heavy states are almost decoupled in low-energy 
effective theory. Alternatively, TeV-scale right-handed neutrinos are also viable, 
which in turn means much smaller couplings to the SM sector and their signals would 
not be captured in future collider experiments such as the Large Hadron Collider 
(LHC). It is therefore difficult to simultaneously realize tiny neutrino masses and 
detectably large interactions among right-handed neutrinos and the SM particles 
within the framework of four-dimensional Type I seesaw. 

In the previous work by a part of the authors [1], it was pointed out that 
the difficulty is overcome by a simple extension of the SM. We have considered 
a five-dimensional theory where all the SM fields are confined in a four-dimensional 
space-time while right-handed neutrinos propagate in the bulk of extra-dimensional 
space [5|,[6]. With bulk Majorana mass, TeV-scale right-handed neutrinos can gener- 
ate a tiny scale of neutrino masses through the seesaw mechanism and simultaneously 
have sizable couplings to the SM leptons and gauge bosons. The previous work fo- 
cused on the tri-lepton signal with large missing transverse energy pp — > t H^i v{v). 
This process is expected to be detectable at the LHC because only a small fraction of 
SM processes contributes to the background against the signals. It was shown that 
the observation of right-handed neutrinos is possible, though it is limited only to 
the lightest Kaluza-Klein (KK) mode of right-handed neutrinos and the mechanism 
itself seems difficult to be confirmed at the LHC. 

In this paper, we will investigate how this scenario can be observed at the Inter- 
national Linear Collider (ILC). The ILC is the future electron-positron linear collider 
for the next generation of high-energy frontier physics. At the ILC, electrons and 
positrons are accelerated by two opposing linear accelerators installed in an about 
30 km long underground tunnel, and are brought into collision with a center of 
mass energy of 500 GeV-1 TeV. The clean experimental environment of the ILC 
due to the electron-positron collider gives us an opportunity to obtain more infor- 
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mation about the property of right-handed neutrinos and their interactions. We 
focus on the two-jet and single lepton signal with large missing energy, e + e~ — > vN 
(N — > £W, W — > qq), which can be efficiently reconstructed at the ILC and allows 
us to reveal not only the lightest but also higher KK-excited states of right-handed 
neutrinos. As a result, the ILC experiment can completely confirm that the scenario 
is based on the extra-dimensional theory. Furthermore, by observing the flavor de- 
pendence of interactions between right-handed neutrinos and SM particles, it is also 
possible to confirm the mechanism to generate the neutrino masses. We performed 
the analysis of the signal in the various mass hierarchies of neutrino masses (the 
normal, inverted, and degenerate cases). 

This article is organized as follows. In the next section, we briefly review the five- 
dimensional theory with right-handed neutrinos. The simulation framework such as 
a representative point in the parameter space of the model and the simulation tools 
used in our study are presented in Section [3J The details of analysis to observe 
the right-handed neutrinos at the ILC are discussed in Section HI where we show 
expected measurement accuracies of the masses of right-handed neutrinos and their 
production cross-sections at the center of mass energy of 500 GeV and 1 TeV. In 
Section we will discuss how the mechanism of neutrino-mass generation can be 
confirmed at the ILC. Section [6] is devoted to summary. 



2 Physics Model 

We consider the five- dimensional theory on S 1 /Z2 with the S^-radius R. The co- 
ordinates are denoted by y) where the y direction is compactified to the line 
segment y = [0, ttR}. The standard-model fields are confined on the four-dimensional 
boundary at y = 0. In addition, we introduce the bulk gauge-singlet fermions 
Mi {i = 1, 2, 3) and assign the even Z 2 parity to their upper components, i.e. 
Ni(—y) = 75-A/i(?/), so that they contain three-generation right-handed neutrinos 
as zero modest The kinetic and mass terms are given by 

C = iJTipNi - l - [a77(M 75 + M')ijAfj + h.c] . (1) 

The conjugated spinor is defined as M c = 7371 J\f . It is easy to also write down the 
Dirac mass rrid if one introduces a Z2-odd function which could originate from some 



: We do not consider the possibility of generation-dependent parity assignment [7]. 
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field expectation value. In this paper, we take = M' = 0, first considered in 
Ref. [3] as a simple example of higher- dimensional seesaw models. We also have the 
mass term between bulk and boundary fields: 

£ m = WiPLTriijLj 8(y) +h.c, (2) 

generated after the electroweak symmetry breaking. The boundary fermions Lj con- 
tain the left-handed neutrinos z/j. Hereafter we take the basis where M is generation 
diagonalized, and will drop other generation indices for notational simplicity. 
The bulk fermions are KK expanded by Majorana fermions ^ R L as 

■A%, y) = J2 x n R {y)P R n^) + E xl(y)P L n(x)- (3) 

n=0 n=l 

The wavefunctions Xr l are normalized so that the KK-mode kinetic terms are canon- 
ical in four dimensions. The low-energy neutrinos (e.g. singlet fermions) come from 
the boundary neutrinos and the KK modes: (u, ty° R , \E^, ^ l L , ty R , $? 2 L , • ■ ■ ) = (v,N). 
By integrating over the fifth dimension, we obtain the Majorana mass matrix in 
four- dimensional effective theory, explicitly given by \{y c N c )PlJ\A(n) +h.c, 
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where the boundary, KK-, and Majorana masses (m„, Mk, and Mrx) are 
m n = x£(0)m, M Rmn = f \x£M X £, 



rR 



ttR 
-R 



M Kmn = dyx^dyxl, M Lmn = dyxTMxl- (5) 

J-ttR J-ttR 

It is noted that M Kmn becomes proportional to 5 mn when Xr l are the eigenf unctions 
of bulk equations of motion, and M Rmn , M Lmn are also proportional to 5 mn for a 
constant (^-independent) mass parameter M due to the normalization conditions. 



2.1 Seesaw and Electroweak Lagrangian 

We further implement the seesaw operation assuming 0{m n ) <C 0(M RLK ) and find 
the induced Majorana mass matrix for three-generations light neutrinos 

My = -M* D M„ X M D . (6) 
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It is useful for later discussions to write down the electroweak Lagrangian in the basis 
where all the mass matrices are generation-diagonalized. The neutrino interactions 
to electroweak gauge bosons are given in terms of these mass eigenstates (u d , N d ): 



C = ±- 
V2 



Wl e^U um P L (u d + VN d ) + h.c. 

+J^T Z ^ + N d V^P L (u d + VN d ), (7) 

A COS u\y 



where W and Z are the electroweak gauge bosons and g is the SU(2) weilk gauge 
coupling constant. The spinors v d are three light neutrinos for which the seesaw- 
induced mass matrix ([6]) is diagonalized 

M v = U* v M d v Ul U v v d = u-M^M^N, (8) 

and N d denote the infinite numbers of neutrino KK modes for which the bulk mass 
matrix Mn is diagonalized both in the generation and KK-mode spaces by a unitary 
matrix Un '■ 

M N = U* N M d N U ] N , U N N d = N + M~ 1 M D u. (9) 

The lepton mixing matrix measured in neutrino oscillation experiments is given 
by L^mns = U\XJ V where U e is the left-handed rotation matrix for diagonalizing the 
charged-lepton Dirac masses. It is interesting to find in (J7J) that the mo del- dependent 
parts of electroweak gauge vertices are governed by a single matrix V defined as 

V = UlM ] D M^*U N . (10) 

When one works in the basis where the charged-lepton sector is flavor diagonal, U u 
is fixed by the neutrino oscillation matrix. 

The neutrinos also have the Yukawa couplings to the Higgs doublet H in the 
four- dimensional boundary, from which the Dirac mass ([2]) is generated; 

C h = fHWP L L6(y)+ h.c, (11) 

where H = eH*. The doublet Higgs H has a non- vanishing expectation value v and 
its fluctuation h(x) in the lower component. After integrating out the fifth dimension 
and diagonalizing mass matrices, we have the Yukawa interaction 

Ch = ^zZ w ^ mnPLU ^ Ud + VN ^ +h - c - ( 12 ) 

n 

and are determined by the mass eigenstates through Eqs. ([8]) and (Q. 
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2.2 Observable Seesaw 



The interactions between heavy neutrinos and SM fields are described by the mixing 
matrix V both in the gauge and Higgs vertices. The 3 x oo matrix V is determined 
by the mass parameters of neutrinos in the original Lagrangian C + C m . The matrix 
elements in V have the experimental upper bounds from electroweak physics. An- 
other important constraint on V comes from the low-energy neutrino experiments, 
namely, the seesaw-induced mass should be on the order of eV scale, which in turn 
specifies the scale of heavy neutrino mass M^. This can be seen from the definition 
of V by rewriting it with the light and heavy neutrino mass eigenvalues 

V = t(M^X(M d N )-^ (13) 

where X is an arbitrary 3 x oo matrix with XX t = 1. Therefore one naively ex- 
pects that, with a fixed order of M d ~ 10 _1 eV and \V\ > 1CT 2 for the discovery 
of experimental signature of heavy neutrinos, their masses should be very light and 
satisfy M$ < keV (this does not necessarily mean that the seesaw operation is not 
justified as M d is fixed). The previous collider studies of TeV-scale right-handed 
neutrinos jH] did not satisfy the seesaw relation ( {TBI and have to rely on some as- 
sumption for suppressing the necessarily-induced (large) mass M v \ for example, the 
neutrino mass matrix must have a singular generation structure, otherwise it leads 
to the decoupling of heavy neutrinos from collider physics. 

A possible scenario for observable heavy neutrinos is to take a specific value of 
bulk Majorana mass jl] so that the lepton number is recovered in low-energy effective 
theory. In this paper we assume that bulk Dirac mass vanishes, but it is easy to 
include it by attaching wavefunction factors in the following formulas. The equations 
of motion without bulk Dirac mass are solved by simple one-dimensional oscillators, 
and the neutrino mass matrices are found 

Tfl 

\ / 2 s "°7iR 

71 

M Kmn = -^Smn , M Lmn = M5 mn . (14) 

From these matrices, we find the seesaw- induced mass matrix M v and the light 
neutrino mixing with heavy modes as follows: 

M u = m 1 - . — — r m, (15) 

2tan(7i\R|M|) v ; 



v = U v i> d H — == m t 



V N* n+1 + V N' 



(16) 
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The heavy mass eigenstate Nj n (N^ n+1 ) are the Majorana fermions with masses 
H — \M\ + \M\). In the seesaw formula, the effect of infinitely many numbers 
of KK neutrinos appears as the factor irR\M\/ tan(7i\R|M|). An interesting case is 
that \M\ takes a specific value \M\ ~ a/R where a is some half integer [5]: the 
seesaw-induced mass M v becomes tiny as a result of KK-mode summation (not 
only suppressed by the Majorana mass scale). On the other hand, the heavy-mode 
interaction is not suppressed unlike the above naive speculation. These facts realize 
the situation that right-handed neutrinos in the seesaw mechanism are observable 
at sizable rates in future collider experiments such as the LHC and ILC. 

In the KK-mode picture, the mass spectrum is vector-like with a half integer a 
and no chiral zero mode exists; for a = 1/2, the mass eigenstates Nd compose the 
Dirac fermions N n {n = 1, 2, • • • ) with masses M n , 

1 1r) — 1 

N n = ^(iVf - M n = JLp. (17) 

As a result, the lepton number is preserved in the KK-mode sector and their contri- 
butions to the seesaw-induced mass vanish. The model given above is an illustrative 
example for accessible seesaw neutrinos. While there are many other possibilities for 
extra-dimensional seesaw, they are supposed to have a common mass matrix struc- 
ture as a key ingredient for tiny neutrino masses and observability, which could be 
seen from the operator analysis in low-energy effective theory. It would therefore 
be reasonable that the above model is used as a representative of extra-dimensional 
neutrinos. With the use of Eq. f ll6p . the weak interaction of KK Dirac fermions, 
which is relevant to the collider study, turns out to be 



2 COS^iy ^— ^ TlRM n \ 5m 

- V] -E-huA P R N n + h.c 

kRv v. 5 M J 



Here 5m = — \M\ characterizes the scale of small neutrino masses. 

Similarly to the seesaw neutrino mass, heavy neutrinos do not give sizable con- 
tributions to lepton-number-violating processes such as the like-sign di-leptons in 
the final states [9]. In the previous work, we have analyzed the LHC signature of 
the above model focusing on the (lepton-number-conserving) tri-lepton signal with 
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Table 1: Three types of neutrino mass hierarchies used in our simulation study: (N), (I), 
and (D) correspond to the normal, inverted, and degenerate mass spectrum, respectively. 



large missing transverse energy jl] (see also [ID]). It was found that the model gives 
enough excessive tri-lepton events beyond the SM background in a wide region of 
parameter space, and the LHC would discover the signs of neutrino mass generation 
and extra dimensions, while the analysis only included the contribution from the 1st 
KK-excited mode. In the following sections, we will perform the ILC study of the 
same setup, in particular, the observation of higher KK neutrino modes and their 
interactions to the SM particles. 



3 Representative Points and Simulation Tools 



3.1 Constraints on Yukawa Couplings 

Before going to discuss representative points used in our simulation study, we sum- 
marize the neutrino mass and mixing matrices which are mandatory to determine 
the flavor structure of Yukawa interaction. The two matrices are parameterized as 
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where s x (c x ) means sin^ x . (cos^). The Dirac and Majorana phases are denoted by 
5 and tpi t 2, respectively. The neutrino mass differences and the generation mixing 
parameters have been measured at the neutrino oscillation experiments [2]. We take 



their typical values; Am 21 
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9xl0^ 3 eV, Am 32 
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eV, S12 — 0.56, s 23 ~ 0.71, and Si 3 < 0.22. The mass spectrum is allowed to have 
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Table 2: The representative points for ?7mns-(N), (I), and (D) correspond to the normal, 
inverted, and degenerate mass spectrum, respectively. 

three types of hierarchies shown in Table [TJ where we define m tot = (0.67 eV — 
2Am 2 i — Am32)/3, considering the cosmological bound J2i m vi — 0-67 eV [IT] . 

Since the scenario we are studying also affects several physical observables such 
as the flavor- changing processes of charged leptons, it is important to take account 
of constraints on the neutrino Yukawa coupling to have proper representative points. 
By integrating out all the heavy KK fermions from the Lagrangian (IPS"]) , we obtain 
the following dimension 6 operator 0^ 6 \ which contributes to the flavor- changing 
neutral current; 

{6) = ^f(LH)fm&L), f = ^m 1 YULs, (20) 

where Y is the 3x3 orthogonal matrix which generally comes in reconstructing high- 
energy quantities from the observable ones [12]. That corresponds to the matrix X 
in (TIB"]) . The coefficient of the operator receives phenomenological constraints as 
shown in Ref. [13], and then each component of the Yukawa couplings is restricted 
by comparing theoretical predictions with experimental data. 

3.2 Representative Points 

We choose the representative points with Mj = M x 1, namely the right-handed 
neutrino masses are degenerate in the flavor space, and also assume that Y is a real 
orthogonal matrix. As a result, the operator is found to be 

° i6) = -4r{LH)U M KsMtul m $(&L). (21) 

To satisfy experimental constraints from this operator and become small e-/z com- 
ponent, we take the lepton mixing matrix £7mns as shown in Table [2] for each case of 
neutrino mass hierarchy. 
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Detector Performance Coverage 

Vertex detector 5 b < 5 © 10/pfi sin 3/2 9 (/an) | cos 0| < 0.93 

Central drift chamber 5p t /Pt < 5 x 10~ 5 (GeV/c) -1 | cos9\ < 0.98 

EM calorimeter cr E /E = 17%/ \[E © 1% | cosfl| < 0.99 

Hadron calorimeter o E /E = 45%/y/E © 2% | cos6»| < 0.99 

Table 3: The detector parameters used in our simulation study. 

Interestingly, the coefficient of depends only on the parameter 5m, which 
turns out to be constrained as follows in each pattern of neutrino mass hierarchy: 

5 M > 3.3 eV for (N), (22) 
hi > 4.4 eV for (I), (23) 
5 M > 24 eV for (D). (24) 

We will set 5m to the most optimistic value, namely these lower bounds, in the follow- 
ing simulation study. The new study [H] may be given the stronger constraints. The 
change of the constraints can be put in the change of 5m bound. It is straightfor- 
ward to extend them to larger values because collider signals such as the production 
cross sections of heavy KK neutrinos are simply proportional to I /5m, though the 
discovery of the signal will be difficult. 

The compactification radius R is not relevant to the constraint from the operator 
. It is however limited by the LEP experiment since the Dirac masses of KK 
neutrinos are given by M n = (2n — 1)/2R. It is easy to confirm that the constraint 
is not so severe if 1/R > 200 GeV, and we thus use the value 

l/R = 300 GeV, (25) 

as a representative point of 1/R in our simulation study. 

3.3 Simulation Tools 

The signal and SM events have been generated by Physsim [15J. The initial-state 
radiation and bremsstrahlung have been included in the event generations. The 
beam energy spread was set to be 0.14% for electron and 0.07% for positron beams. 
The finite crossing angle between the electron and positron beams was assumed to 
be 14 mrad. In the event generations, helicity amplitudes were calculated using the 
HELAS (HELicity Amplitude Subroutines) library [16] , which allows us to deal with 
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Figure 1: The diagrams of signal processes, e + e — > uN (N — > IW). 

the effect of gauge boson polarizations properly. The phase space integration and the 
generation of parton 4-momenta have been performed by a nubmerical integration 
program package , BASES/SPRING [17]. The parton showering and hadronization 
have been carried out by using the event generator PYTHIA6.4 [18] . where final-state 
tau leptons are decayed by the MC particle decay software package, TAUOLA [19] 
to handle their polarizations correctly. 

The generated Monte Carlo events have been passed to a detector simulator called 
JSFQuickSimulator, which implements the geometry and other detector-performance 
related parameters of an ILC detector concept called GLD [20] . In the detector sim- 
ulator, hits by charged particles at the vertex detector and track parameters at the 
central tracker are smeared according to their position resolutions, taking into ac- 
count correlations due to off-diagonal elements in the error matrix. Since calorimeter 
signals are simulated in individual segments, a realistic simulation of cluster over- 
lapping is possible. Track-cluster matching is performed for the hit clusters in the 
calorimeter in order to form pseudo Particle Flow Objects (pPFO), thereby achieving 
the best attainable jet energy measurements. The resultant detector performance in 
our simulation study is summarized in Table [3j 

4 Results from Simulation Study 

The simulation has been performed at y/s = 500 GeV for the first KK mode of 
right-handed neutrinos, and at a/s = 1 TeV for the KK modes up to the third level, 
with an integrated luminosity of 500 fb _1 each. We have considered the case with no 
beam polarization. We have evaluated the measurement accuracies of the masses of 
KK right-handed neutrinos and their production cross-sections using the processes 
as shown in Fig. [1] for the cases of the normal, inverted, and degenerate hierarchies 
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Table 4: Cross sections of e+e" -> z^iV (iV -»• £W, W -> gg) at ^/i = 500 GeV and 1 TeV. 
The three dots mean that the cross sections are too low to be explored at the ILC. 

of neutrino masses. 



4.1 Study at = 500 GeV 

Based on the cross sections shown in Table HI we have studied e + e~ — > uN% (N± — > 
eW) for all the neutrino mass hierarchies and e + e~ — > uN\ (N\ — > rW) for the 
normal and inverted mass hierarchies. In the analysis, we have used the hadronic- 
decay modes of W, which allow us to fully reconstruct the mass of Ni. 

4.1.1 Analysis of vNi -» veW 

In the signal event, an isolated-electron track from the decay of N\ is expected. We 
have therefore selected the electron track and reconstructed two jets in the final state 
of the signal. Since the isolated-electron track has no energy around it, while a track 
from a jet has some energy, we have selected tracks with the energy around them 
within 20 degrees below 5 GeV. Then, in the tracks which satisfy the requirement, 
the track with the maximum energy has been selected as a candidate of the electron 
track. After picking up the isolated-electron track, the clustering of the jets has been 
performed. The pPFOs have been combined to form a jet if the two clusters satisfy 
Uij < Vcut, where the variable is defined as 

2E i E i (l-cos6 ii ) . , 

Vij = A ^. (26) 

Here, 6^- is the angle between two clusters, are their energies, and E vis is the 
total visible energy. All events are forced to have two jets by adjusting y cut . The 
mass of Ni was, then, reconstructed using the candidate for the electron track and 
two reconstructed jets. 
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Figure 2: Distribution of the N± mass reconstructed from M e jj in the case of inverted 
neutrino mass hierarchy at y/s = 500 GeV. 

We have considered the processes shown in Table [5] as the background. In order 
to suppress these background processes, we have applied the following requirement. 
The isolated electron candidate selected above would have a higher energy then 
tracks from jets, if it were really coming from the decay of Ni. We, hence, required 
the energy of the isolated-electron candidate (E e ) to be 10 GeV < E e < 200 GeV. 
In addition, since the reconstructed di-jet mass {Mjj) should be consistent with the 
W hypothesis for a signal event, we have selected events with 60 GeV < Mjj < 
100 GeV. For reconstruction of the Ni mass with the isolated electron and W 
candidates (M e jj), the energy and momentum of the W candidate was corrected to 
have W mass. Figure H] shows the distribution of M e jj after all the selection cuts in 
the case of the inverted neutrino mass hierarchy. With the signal region defined by 
135 GeV < M e jj < 165 GeV, the numbers of signal and background events before 
and after selection cuts are summarized in Table [51 

We have evaluated the measurement accuracy of the Ni mass and its production 
cross section by fitting the M e jj distribution. The Ni mass values and errors from 
the fit are 150.2 ± 0.20 GeV, 150.0 ± 0.01 GeV, and 150.0 ± 0.01 GeV, for the 
normal, inverted, and degenerate neutrino mass hierarchies, respectively. These 
results correspond to the measurement accuracies of 0.14%, 0.01%, and 0.01%. On 
the other hand, the cross section of e + e~ — > vNi [N\ — > eW) events turned out to 
be measurable with accuracies of 6.5%, 0.4%, and 0.4%. Since the cross-section of 
the normal hierarchy is much smaller than that of the other hierarchies as shown in 
Table HI the measurement accuracy for the normal hierarchy is worse. The results 
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are summarized in Table [7J 
4.1.2 Analysis of z/Aq — > vtW 

In order to analyze e + e~ — > vN\ (Aq — > tW) events, we have reconstructed all events 
as 3-jet. A jet with the smallest number of the tracks was assumed to be a tau-jet. 
The previous signal e + e~ — > z/Aq (Aq — > eW) has been considered as background 
together with e + e~ — > evW and WW shown in Table [5j The following selection cuts 
have been applied to reduce these backgrounds. The energy of the tau-jet (E T ) was 
required to satisfy 10 GeV < E T < 150 GeV in order to reject high energy electron- 
and muon-tracks from the leptonic-decay modes of W in the e + e~ — > WW events. 
We required the same criteria for the di-jet mass, Mjj, as in the study of uN — > eW. 
Then, the reconstructed Aq mass (M T jj) was required to be 80 GeV < M T jj < 160 
GeV. 

After applying the selection cuts, a likelihood analysis has been performed. Since 
the processes e + e~ — > euW, z/Aq (Aq — > eW) dominate in the background, we 
constructed two likelihood functions to separate the signal from evW (£ e vw) and 
z/Aq — > veqq {C- V N x ^veqq)- As the input variables of the likelihood functions, we used 
the number of tracks in the jets of the tau candidate, the energy of the track with 
the maximum energy in the tau-jet (£ maj! ), and the energy of tau-jet with ^ 
subtracted from it. These likelihood functions were prepared for the normal and 
inverted neutrino mass hierarchies separately, because we would be able to identify 
the neutrino mass hierarchies by using the cross section of uN — > veW events. 
We have required C eu w > 0.79 (0.63) and C^N^ueqq > 0.13 (0.11) for the normal 
(inverted) neutrino mass hierarchy to maximize the signal significance. 

The resolution of the Aq mass can be improved by compensating for the missing 
energy of the r decay as follows: Since the Aq mass has already been measured 
through the analysis of uN — > veW events, we can calculate energy of the Aq 
assuming its two-body kinematics. Then, we calculated the energy of the tau as the 
calculated Aq energy minus the energy of the W candidate. We assumed that the 
direction of the tau coincided the direction of the tau-jet. The corrected mass of Aq 
( col M T jj) was reconstructed by using the estimated tau energy and momentum and 
those of the W candidate. Figure |3] shows the distribution of col M T jj for the inverted 
neutrino mass hierarchy after applying all selection cuts. The number of events in 
the signal region, 135 GeV < co1 M T jj < 165 GeV, before and after the selection cuts 
are summarized in Table |5j Fitting the col M T jj distribution, we have obtained the Aq 
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Figure 3: Distribution of the N± mass reconstructed from co[ M T jj in the case of inverted 
neutrino mass hierarchy at \fs = 500 GeV. 

mass as 149.8 ± 0.24 GeV and 150.0 ± 0.32 GeV, corresponding to the measurement 
accuracies of 0.16% and 0.21% for the normal and inverted neutrino mass hierarchies, 
respectively. The cross section of vN\ — > vrW could be determined with accuracies 
of 11.3% and 12.4% for the two mass hierarchies. These results are summarized in 
Table [3 

4.2 Study at = 1 TeV 

Not only the first but also the second and the third KK modes of right-handed 
neutrinos can be produced at the ILC with ^fs = 1 TeV. Taking into account 
the cross sections shown in Table HI we have studied the process e + e~ — > vNi 
(Ni — > veW) for all the neutrino mass hierarchies and e + e~ — > vN 2 ^ (^2,3 — > veW) 
processes for the inverted and degenerate neutrino mass hierarchies. The background 
processes considered in this study are shown in Table EJ 

After the selection of the electron track and the reconstruction of two jets with 
the same procedure as the study at y/s = 500 GeV, the following selection cuts were 
applied. To remove electron tracks from jets, we selected high-energy electrons by 
requiring 10 GeV < E e < 600 GeV. The di-jet mass was required to satisfy the same 
criteria as the study at yfs = 500 GeV. Since the WW events have a peak at 500 
GeV for the di-jet energy distribution, we have required the di-jet energy should be 
below 460 GeV. Shifting the energy and momentum of the W candidate to have 
the nominal W mass, the M e jj was calculated. Finally, the signal region has been 
defined as 135 GeV < M ejj < 165 GeV, 425 GeV < M ejj < 475 GeV, and 720 GeV 



I I vNi ->-vtW 
vNi-»veW 
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Figure 4: Distributions of the masses of (a) the second and (b) third KK modes recon- 
structed from M e jj in the case of inverted neutrino mass hierarchy at y/s = 1 TeV. 

< M e jj < 780 GeV for the 1st, 2nd, and 3rd KK modes, respectively. The numbers 
of events before and after the selection cuts are summarized in Table [6j 

Figure H] shows the distributions of M e jj for the 2nd and the 3rd KK modes of 
right-handed neutrinos in the case of inverted neutrino mass hierarchy. Fitting the 
M e jj distributions, we have obtained the measurement accuracies of iVi^^ masses 
and their cross sections as shown in Table The masses of right-handed neutrinos 
could be determined with the accuracy better than 1%. 

5 Discussion 

We found that the observation of the higher KK modes, especially the mass spectrum 
of these particles is measured accurately at the ILC, which allows us to confirm that 
physics behind these signals is based on higher- dimensional theory. The observation 
of the masses of KK modes, however, does not directly mean that the signals are 
coming from physics responsible for neutrino masses and mixings, because there are 
many scenarios in the framework of higher-dimensional theory which are not related 
to neutrinos but to other issues such as hierarchy problem and dark matter. 

In order to confirm that the signals are from physics of neutrino masses and mix- 
ings through the higher- dimensional theory, we should observe not only the masses 
of KK modes but also other quantities which depend strongly on parameters of 
neutrinos. One of such parameters is 5m, because it is the origin of tiny neutrino 
masses. However, the smallness of the parameter inevitably leads to the smallness of 
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Process (</s = 500 GeV) 


Cross sec. 


Events No. 


Events after cuts 


uNi ->■ veW (W ->■ gg) (N) 

(I) 


6.5 [fb] 
297.5 [fb] 

ZOrf.l [IDJ 


3,262 
177,700 

lZo,0D4 


2,216 
101,057 

87 coo 
o 1 ,Ooo 


ez/IV — >■ ez/gg 
WTV -> £i/gg 

— > z/i/gg, ££gg 


4,462 [fb] 
1320 [fb] 
108 [fb] 
531 [fb] 


2,231,000 
660,000 
54,000 
265,500 


22,594 
11,324 
1 

568 




Process (^/s = 500 GeV) 


Cross sec. 


Events 


Events after cuts 


i/iVi -»■ vtW (W -+qq) (N) 

(I) 


5.49 [fb] 
4.18 [fb] 


2,745 
2,090 


1,029 
821 


vN± -»■ rely (IV -> gg) 
ez/IV — >■ ez/gg 
WW -> £z/gg 


6.52 [fb] 
4,460 [fb] 
3,960 [fb] 


3,260 
223,100 
1,980,000 


495 (N) 554 (I) 
8,989 (N) 12,276 (I) 
13,788 (N) 14,861 (I) 



Table 5: Summary of cuts for y/s = 500 GeV. 



lepton-number violation in the scenario, so that it is difficult to observe the quantities 
related to this parameter through processes violating the lepton-number. 

Another interesting parameter is £/mns describing phenomena of neutrino mix- 
ings. The higher-dimensional theory for neutrinos considered in this article predicts 
a specific flavor-structure at the right-handed neutrino sector, and it is completely 
determined by the mixing parameter C/mns an d neutrino masses Mf, as shown in 
Eq. ([TBI) . This fact leads to that branching fractions of the right-handed neutrino 
are governed by Umns an d M$. We have shown that both the eW and rW de- 
cay modes of Ni can fortunately be measured at the ILC. The ratio of these two 
branching fractions R = Br(rW)/Br(eTV) is theoretically given by the formula, 

i2 r , 1 r 

^MNS {M^) 2 U MNS 



R 



U- 



MNS 



LVL v ^MNS 



13 



u 



MNS 



33 



11 



MNS 



M d ift 

lvl u U MNS 



£/mns {M$y ui 



MNS 



(27) 



33 



which takes a value of 0.850 for the normal mass hierarchy and 1.58 x 10 -2 for the 
inverted mass hierarchy using the representative points of C/mns m Table [2J With the 
use of the results in Table H, we investigate how accurately the ratio can be measured 
at the ILC. The result is shown in Fig. [5] for both cases of normal and inverted mass 
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Process (y/s = 1 TeV) 


Cross sec. 


Events 


Events after cuts 


i/JVi ->■ i/elf ->■ qq) (N) 

(I) 
(D) 


7.79 [fb] 
355.0 [fb] 
307.0 [fb] 


3,895 
177,700 
153,540 


1,244 
44,021 
29,106 


ez/VT — > euqq 
WW ->■ £z/gg 

— >■ z/z/gg, ^£gg 

ft 


10,320 [fb] 
560.6 [fb] 
42.79 [fb] 
29.43 [fb] 


5,160,000 
280,300 
21,393 
14,715 


9238 
1,234 
5 
23 




Process (^s = 1 TeV) 


Cross sec. 


Events 


Events after cuts 


vN 2 ->■ z/eW -> gg) (I) 

(D) 


23.6 [fb] 
20.4 [fb] 


11,800 
10,200 


6,756 
5,820 


ez/VT — >■ ez/gg 
WW -> £z/gg 
ZZ — >■ z/z/gg, ^£gg 


10,320 [fb] 
560.6 [fb] 
42.79 [fb] 
29.43 [fb] 


5,160,000 
140,150 
21,393 
14,715 


24,671 
4,858 






Process (-^/i = 1 TeV) 


Cross sec. 


Events 


Events after cuts 


z/iV 3 ->■ z/eW ->■ qq) (I) 

(D) 


3.86 [fb] 
3.34 [fb] 


1,932 
1,670 


1,131 
961 


ez/VT — > euqq 
WW ->■ £z/gg 

— >■ z/z/gg, ££gg 

ft 


10,320 [fb] 
560.6 [fb] 
42.79 [fb] 
29.43 [fb] 


5,160,000 
280,300 
21,393 
14,715 


10,510 
6,780 





Table 6: Summary of cuts for yfs = 1 TeV. 



17 



Process (y/s = 500 GeV) 


iV mass resolution 


Cross-section accuracy 


vNi ->■ veW (N) 

(I) 
(D) 


0.14 [%] 
0.01 [%] 
0.01 [%] 


6.5 [%] 
0.4 [%] 
0.4 [%] 


vN x ->■ i/rlf (N) 

(I) 


0.16 [%] 
0.21 [%] 


11.3 [%) 

12.4 [%] 




Process (a/s = 1 TeV) 


mass resolution 


Cross-section accuracy 


i/JVi ->■ rely (N) 

(I) 
(D) 


0.41 [%] 
0.01 [%] 
0.02 [%] 


13.6 [%] 
0.6 [%] 
0.7 [%] 


vN 2 ->■ i/eiy (I) 

(D) 


0.05 [%] 
0.08 [%] 


2.8 [%] 
3.1 [%] 


z/iV 3 ->■ i/eiy (I) 

(D) 


0.21 [%] 
0.23 [%] 


9.9 [%] 
10.0 [%] 



Table 7: Summary of measurement accuracies. 




R R(x10") 



Figure 5: Probability Densities to determine the branching ratio between the eW and tW 
decay modes of the right-handed neutrino N±. The left figure shows the density for the 
case of normal mass hierarchy, while right one is for inverted mass hierarchy. 
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hierarchies, where probability densities for the determination of the ratio are shown. 
It can be seen that the ratio will be measured with accuracy of ± 0.110 and ± 0.17 
x 10~ 2 for the cases of normal and inverted mass hierarchies, respectively. The ILC 
has, therefore, a potential not only to discover the signal of higher- dimension theory 
but also to observe the signal related to neutrino oscillation parameters. 

6 Summary 

We have investigated ILC signals of the seesaw scenario in a five-dimensional exten- 
sion of the SM, where right-handed neutrinos live in the bulk and the SM particles 
stay at a four-dimensional boundary. We focused on the production process of KK 
right-handed neutrinos, e + e~ — > Nv (N — > £W, W — > qq), where the masses of 
KK right-handed neutrinos N can be fully reconstructed. With realistic Monte- 
Carlo simulations, we found that the masses of KK neutrinos and their production 
cross sections can be measured accurately at the ILC as summarized in Table [7J In 
particular, it was found that the mass and production cross section of the first KK 
right-handed neutrino can be measured accurately for various hierarchies of neutrino 
masses with the center of mass energy of 500 GeV. In addition, it was shown that 
masses and production cross sections of the second and third KK neutrinos can be 
measured with the center of mass energy of 1 TeV. 

The ILC also allow us to investigate the flavor structure of the higher-dimensional 
theory using the measurement of the N± production cross section followed by its 
various decay modes. Since branching fractions of N\ decay modes are determined by 
the masses of neutrinos and their mixing matrix C/mns; measuring the ratio between 
these fractions can be directly compared with the results of neutrino oscillation 
experiments, which give us an important clue to clarify the mechanism to generate 
neutrino masses in the framework of the higher dimensional theory. Interestingly, 
the first or third generation lepton in the final state of the Ni decay can be detected 
at the ILC, while the LHC will be possible to observe that with the first or second 
lepton. The ILC will therefore be a complementary machine to the LHC to explore 
physics of neutrinos in the higher- dimensional theory. 
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